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Abstract 

We extend a recent study that explained the steep falloff in the fusion cross section at energies far 
below the Coulomb barrier for the symmetric dinuclear system ^^Ni+^^Ni to another symmetric 
system, ^^Ni+^'^Ni, and the asymmetric system ^^Ni+^'^'^Mo. In this scheme the very sensitive 
dependence of the internal part of the nuclear potential on the nuclear equation of state determines 
a reduction of the classically allowed region for overlapping configurations and consequently a 
decrease in the fusion cross sections at bombarding energies far below the barrier. Within the 
coupled-channels method, including couplings to the low-lying 2"*" and 3~ states in both target and 
projectile as well as mutual and two-phonon excitations of these states, we calculate and compare 
with the experimental fusion cross sections, S'-factors, and logarithmic derivatives for the above 
mentioned systems and find good agreement with the data even at the lowest energies. We predict, 
in particular, a distinct double peaking in the 5-factor for the far subbarrier fusion of ^^Ni+^^Ni 
which should be tested experimentally. 
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I. INTRODUCTION 



Very recently we proposed a mechanism that can explain the unexpected hindrance or 
steep falloff of fusion cross sections which has been observed at bombarding energies far 
below the Coulomb barrier Although measurements of several medium-heavy nuclei 
performed in the past two decades already provided some indications of a steep decrease 
of the excitation functions at the lowest bombarding energies, the credit of disclosing and 
confirming this unexpected trend for new fusing systems has to be given to C. L. Jiang et 
al. [2, 3], y, [5 , 0] . Among the most conspicuous cases reported in the past are ^^Ni+^^Ni [3], 
where the departure from the expected behavior takes places already at cross sections ~ 0.1 
mb, and ^°Zr+^^Y, ^°Zr+^^Zr as reported in j^. 

The new fusion data reported by Jiang et al. are even more spectacular because the 
reported cross sections are measured down to 10 nb : ^°Ni+ ^^Y 121 (af > 100 nb), ^^Ni-|- 
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6^Ni y {aj > 10 nb), 6^Ni+i°°Mo (sj (cx/ > 10 nb). The hindrance of fusion was first 
reported as a suppression of the measured low-energy fusion cross sections with respect to 
model calculations It was later on characterized by the energy Es where the S-factor for 
fusion develops a maximum at low energy . Following the publication of these findings a 
challenge was launched, especially on the theoretical side of the heavy-ion fusion community. 
At the end of 2005 the underlying physical cause of this apparently new phenomenon was 
still unknown according to the authors of Ref. |[6j. Some authors have even advocated 
the hypothesis that the standard theoretical approach to treat capture reactions, i. e. the 
coupled-channels (C. C.) method is unable to explain the steep falloff of the cross sections. 
We shall try to convince the reader throughout this paper that the C. C. method is the right 
tool to investigate capture, even at very low energy. 

The nuclear potentials that are employed in C. C. calculations are commonly parametrized 
as a Woods-Saxon well. Among the issues related to the deep sub-barrier fusion was the 
large diffuseness a of the ion-ion potential that was needed to fit high-precision fusion data. 
Hagino et al. [ol, hinted that a phenomenological nuclear potential with a larger diffuseness 
leads to a better agreement to the data. Values up to a = 1.3 fm were conjectured instead 
of the usual a = 0.65 fm for the system ^^Ni+^^Ni. In it was remarked that since 
the low-energy fusion becomes sensitive to the nuclear potential inside the barrier, this 
part of the interaction may not be accurately modeled by the conventional Woods-Saxon 



2 



parametrization. It was pointed out that by doubling the diffuseness of the inner part of 
the potential, the agreement with the data improves for the colliding system ^''Ni+^^Y sl- 

Following the same idea, the systematic failure of the Woods-Saxon nuclear potential to 
describe fusion was analyzed in 10| and it was concluded that the origin of it should be 
sought in the diffuseness parameter a. In order to fit the data at energies above the Coulomb 
barrier, the diffuseness must be increased with increasing Z1Z2. The fusion data compiled 
by these authors indicate a correlation with the neutron richness of the projectile and target 
nuclei in the sense that the neutron rich nuclei tend to require larger values of a. 

The authors of [HI] pointed out that potentials such as the Akyiiz-Winther (AW) 
are providing reliable barriers but they cannot reproduce the data far below the barrier, a 
fact which made them suggest that the ion-ion potential should have another form in the 
inner part of the barrier. Following a sequence of simple but clear arguments they pointed 
out that the exponential falloff in the tunneling probability is related to the disappearance 
of the classically allowed region below a certain energy. If this is true, it would mean that 
we are confronted with the existence of a shallow pocket of the potential inside the barrier. 

In Ref. IJ] a surprisingly good description of the data for ^^Ni+^^Ni, ^"^Ni-l-^^Ni and 
^°Ni+^^Y was claimed. However, it is difficult to judge the significance of the results be- 
cause of the limited number of excitations that were included in the C. C. calculations. For 
example, the only excitation considered in the nickel isotopes is the one-phonon 2+ exci- 
tation, but it is well known that couplings to the 3~ state and higher-order cou plin gs to 
two-phonon states are tremendously enhancing the heavy-ion fusion cross section [15|, 

Other works have attempted to account for the complicated problem of channel coupling 
by means of a polarization potential 17|. The imaginary potential parametrizes in this case 
;he excitation of other degrees of freedom that influence the fusion process. According to 



16|. 



171] , the imaginary potential shows a rapid cutoff at energies far below the Coulomb barrier 



for the cases ^^Ni+ ^^Ni, ^^Ni+ '^^Ni and '^^Ni+ ^"^Ni because of a threshold below which the 
C. C. effects are ceasing to exist. 

Before ending the list of hypotheses advanced to explain the fusion hindrance phenomenon 
we quote the result reported in Ref. jisl, which overrules the possibility of explaining the 
depletion of fusion rates at extreme sub-barrier energies from dissipative tunneling, as results 
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from other quantum open system approaches, such as that of Caldeira and Leggett 

Very recently |l| we proposed an explanation of the hindrance observed in the sub-barrier 
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fusion of 6^Ni+^^Ni, which was based on the same standard coupled- channels formalism as 

n 

before p| but with amendments that concern the potential. Essential in getting a good 
description of the data was to take into account the saturation of nuclear matter and to use 
realistic neutron and proton distributions of the reacting nuclei. These two ingredients are 
naturally incorporated in a potential calculated via the double-folding method with tested 
effective nucleon-nucleon forces and with realistic charge and nuclear densities, a fact which 
is often overlooked or only indirectly included in the Woods-Saxon parametrization. We 
intend in this paper to show that by properly addressing these issues, light can be shed on 
the extreme sub-barrier fusion data. Before we do that we summarize the methods that 
have been used to analyze the low-energy behavior of heavy-ion fusion cross sections. 



II. REPRESENTATIONS OF LOW-ENERGY CROSS SECTIONS 



In an attempt to get a diagnosis for the various cases where the hindrance in sub-barrier 



fusion occurs, the authors of Ref. [3| proposed the use of two representations. The first one 
is the astrophysical S'-factor [2^, 

S = Ect{E) exp{2nr]), (1) 

where E is the center-of-mass energy, rj = ZiZ2e^ / {hv-tci) is the Sommerfeld parameter and 
frei is the relative velocity of the fragments. The experimental value of S increases with 
decreasing bombarding energy and has the tendency to develop a maximum for the systems 
of interest. The necessity to resort to this quantity comes from the fact that the reaction 
cross section varies by many orders of magnitude below the Coulomb barrier (7 orders of 
magnitude for ^^Ni+ ^"^Ni). The S'-factor has been used in the past to unravel typical 
molecular resonant structures in the excitation function of systems like ^^C+^^C since it 
removes the dominating influence of the Coulomb and centrifugal barrier transmission factors 



that mask these structures in the cross section 2l|. The series of narrow and prominent 
resonances was associated with quasi-bound, long-lived states of the ^^Mg nucleus. 

Very recently the fusion cross section for the ^^C+^^C system has been measured down 



to very low energies 



22| . The data show a rise in the S'-factor at the lowest energies, which 



might indicate the existence of a broad resonance in the entrance channel, possibly related 



to an intermediate state in the compound nucleus 



22| . Similar broad resonances in the 
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S'-factor at energies below the barrier have been also inferred from the C+ O total cross 



sections 23|]. Thus, the S'-factor is a quantity that magnifies structures in the excitation 
function at energies below the barrier, and it is also an instrument for exploring the inner 
part of the barrier in low-energy, heavy-ion fusion reactions. 

A second representation proposed in Ref. is the logarithmic derivative. 

The point where the experimental L{E) intersects the logarithmic derivative obtained from 
an s-wave transmission across a pure point-charge Coulomb potential (constant S'-factor), 
given by Lq{E) = coincides with the maximum in the S-factor invoked earlier. Extracting 
the energy Eg, where this intersection occurs, it was found that the corresponding logarithmic 
derivative Lg is nearly the same for stiff heavy-ion systems, with an average value of 2.34 
MeV~^. This implies that Es scales with the charge and mass numbers of the reacting nuclei 
according to the empirical law (sl 

2/3 



E, = 0.355 



Z1Z2 



A1A2 



(3) 



Ai + A2 

At this point it is useful to recall that the inability of previous calculations to repro- 
duce the low-energy data points of the measured fusion cross sections is most clearly seen 
from the inspection of the logarithmic derivative L{E). Thus for energies below a certain 
threshold, the experimental values of L{E) increase steeply with decreasing energy, whereas 
the theoretical curve increases with a much smaller slope on which a resonant structure is 
superposed. 

III. SIMULATIONS OF A REPULSIVE CORE 

The occurrence of resonant structures in collisions of light nuclei, the best known example 
being the sharp peaks in the bombarding energy dependence of the gamma radiation yields 



emerging in the ^^C+^^C scattering found by Bromley et al. 2J], were found to resemble 
states in a molecular potential well. Following the suggestion that these " quasi- molecular" 
states may represent doorway states to fusion, the resonant behavior found its explanation 
by introducing the concept of the double-resonance mechanism [25]. In this scenario the 
indirect population of quasi-molecular states can occur in light-ion scattering according to 



the following sequence of events: a) Surpassing the potential barrier at an initial energy Ei 
and losing the energy E* by inelastic excitations of low-energy levels of one or both of the 
ions, and b) Resonant decay into the potential pocket, where the colliding nuclei are trapped 
and a quasi-molecule is formed, if their relative energy E.^ — E* coincides with the energy of 
a quasi-bound state. 

Evidence that the resonant behavior observed in ^^C+^^C is not an isolated phenomenon 
was made available for other light systems such as ^^C+ ^^O, ^^C+ ^^C, ^^0+ ^^O, ^^0-|- 
^^Mg. It persists even in heavier colliding systems, such as the ^'*Mg-|- ^^Mg 26|] and ^^Si-|- 
^^Si [27]. Manifestation of clusterization in connection to quasi-molecular pockets is also 



known for heavy nuclear systems, such as the cluster radioactivity 28| , or hyperdeformation 



and clustering in the actinide region 
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To simulate the appearance of shallow pockets several recipes have been proposed : 

• A central soft repulsive core added to the conventional Woods-Saxon potential was 
used in Refs. jsi] and 32] in order to fit the reaction cross sections observed in ^^C-|- 
i2q^ i2q_|_ i6q 1^0+ ^^O. The effect of a repulsive core was modeled by a positive 
Gaussian potential, Kep exp(— 6r^), the width h being constrained by the requirement 
that the potential becomes repulsive for r < Rq, Rq being the radius of the Woods- 
Saxon potential. 



Double-folding potentials as introduced by Satchler and Love 33|] are accurate only in 
the tail region of the nucleus-nucleus potential, where the density distributions are only 
gently overlapping and thus the assumption of "frozen density" is less questionable. 
However this assumption ignores any readjustment due to the mutual interaction of the 
nuclei or the Pauli exclusion principle for strong overlap. To cope with this problem 



the observation was made in 



3^, l35(] that, whereas any theory of heavy-ion potentials is 



expected to reduce to the double-folding model in the limit of large ion-ion separations 
and vanishing density overlap, the compound system resulting from fusing the two ions 
is accurately described by the liquid-drop energy -Eldm- To interpolate between these 
two extremes the nuclear deformation energy was written as 

Edef = -ElDM + Wm3Y+c{p) - (yWM3Y{p), (4) 

where Wm3y+c{p) is the double folding self-interaction energy with the original M3Y 
parametrization for the nuclear and Coulomb forces with realistic values of the proton 



and neutron matter diffusivities and radii. The Wm3y{p) is a renormalization intro- 
duced with the purpose to subtract the volume and surface energy contributions in 
W, and is computed with sharp-edge densities. The parameter a is adjusted so that 
W and W cancel exactly in the limit of complete overlap. This scheme was applied 
only for collisions of two identical spherical nuclei since the evaluation becomes to 
cumbersome for deformed nuclei approaching at different orientations. 

Pockets in a double folding potential are arising naturally if an effective Skyrme like 
A^A^ force is used [36|]. However, this force was not tested systematically for scattering 
reactions as is the case for the M3Y force. 

Density-dependent interactions superposed on the original M3Y form have been used 
in recent times to simulate the saturation of nuclear matter for a+nucleus scattering 



or elastic scattering of light nuclei |37|, 



381]. An attempt to explain fusion data for a 



variety of systems using these density dependent interactions was made very recently 



by the Canberra group 39|] and they reached the conclusion that the double-folding 
model may not be appropriate for fusion. The pockets resulting from the density 
dependence of the effective A^A^ force are still too deep to improve the agreement of 
C. C. calculations and the experimental data at extreme sub-barrier energies. 

Methods based on the energy density functional are known for a long time to be 

n n 

able to predict shallow pockets in the interaction between two nuclei [40, |41|. A 
typical feature of this class of potentials is the short-ranged repulsion at distances 
where a strong overlap of the nuclear densities takes place. In this approach the 
condition of nuclear matter saturation is achieved by fitting the free parameters in the 
energy- density functional using the properties of finite nuclei. Like in the case of the 
double-folding recipe, the sudden-approximation, i.e. the summation of frozen nuclear 
densities, is essential for the occurrence of the core. Very recently, using this method, in 
the framework of the extended Thomas-Fermi approximation with fi^ correction terms 
in the kinetic energy density functional, an analytic potential with parameters that 



are fitted to data was proposed 42] . However, this potential provides barriers that are 
higher than the experimental ones and they can therefore not provide a satisfactory 
description of the data in the barrier region. 
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• Proximity potentials are well-known examples of nuclear interactions producing pock- 
ets in the ion-ion potential (see Ref. j43| for the 1977 and [44] for the 2000 versions). 
Unfortunately these potentials provide pockets that are too deep for the systems of 
interests (see for example Fig. 1 of Ref. 

• Finally we mention the semi-empirical Akyiiz-Winther potential which we will refer to 
in the following as the AW potential. It is parametrized as a Woods-Saxon potential 

12|, so that the maximum nuclear force is consistent with the proximity force for 
touching spheres, and it has an exponential tail that is consistent with the M3Y 
double-folding calculated in Ref. The AW also produces pockets that are too 

deep as evidenced by its applications in Refs. [l, 2, ^, 0, Q]. 

In the present work we are interested in the interaction between two intermediate mass 
nuclei. For these target-projectile combinations the above approach to the density dependent 
interaction does not necessarily apply because, as we shall see, the resulting pockets are still 
to deep. For this reason another approach will be adopted. 



IV. THE M3Y-KREPULSION, DOUBLE-FOLDING POTENTIAL 

As a manifestation of the Pauli principle, which prevents the overlapping of the wave 
functions of two systems of fermions, we expect that the interaction potential between two 
colliding nuclei will contain a repulsive core. The actual form of the repulsive core and 
its strength depend strongly on the extent to which the collision is adiabatic or sudden. 
A further uncertainty with regard to the core parameters is the influence of individual 
characteristics of the considered nuclei, including binding energies, shapes, and the nucleon 
distributions. The models of ion-ion potentials which provide a repulsive core and which we 
enumerated in the preceding section, lead to quantitatively different estimates of the height, 
radius and diffuseness of the core potential. Due to these conditions of strong uncertainty 
we use like in previous works on fusion a crude recipe to determine the properties of the 
repulsive core potential. 
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A. Calibration of the repulsive core potential 



As noticed earlier, an overlapping region with doubled nucleon density is formed once the 
distance r between the nuclei becomes less than Ri + R2, where Ri and R2 are the nuclear 
radii of the target and projectile along the scattering (fission) axis. Here and in the following 
we assume for simplicity that the densities of the two nuclei are frozen when estimating the 
nuclear interaction potential. When a complete overlap of the two nuclei occurs, the total 
density is therefore roughly twice that of normal matter, p ^ 2po, within the volume of the 
smaller nucleus. Consequently, the nuclear equation-of-state (EOS) dictates an increase AV 
in the energy of the compound system. In the case of complete overlap (for r = 0) AV is 
proportional to the increase of the energy per particle of nuclear matter e{p, 6) considered as 
a function of the nuclear density, p = pn+Pp, and the relative neutron excess, 6 = (pn—Pp) / p, 

AV ^2Ap[e{2po,S)-e{po,S)]. (5) 

Here the relative neutron excess 6 is assumed to be approximately constant, and Ap is the 
mass number of the smaller nucleus in the case of an asymmetric system. 



The EOS predicted by the Thomas-Fermi model for cold nuclear matter is 

2/3 / ^ \ / ^ ^ 5/3" 



e{p,6) = Ef 



A{6) 



Po) \PoJ \poJ 



(6) 



where po = 0.161 fm ^ is the saturation density and ep is the Fermi energy of normal nuclear 



matter. The expressions for the (5-dependence entering Eq. ([6]) are listed in Ref. [46||. The 
definition of the incompressibility of cold nuclear matter is then related to the curvature of 
e{p, 6) at the saturation density, 

A-9(,=0) . (7) 

In order to calibrate the strength of the repulsive core potential, we assume that in 
Eq. (IS]) must be identified with the nuclear part of the heavy- ion potential, VAr(r), evaluated 
at the coordinate origin r = 0, 

AV = V^^(O). (8) 



Approximating the EOS by a parabolic expansion around the equilibrium density po |47|, 

K 

e{p,6) =e{po,S) + —^ {p- PoY, (9) 



we obtain the following estimate of the nuclear potential at r = 

A 



VW(0) = ^ ^ K (10) 
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They are 



Similar recipes for introducing a repulsive core can be found in the literature 
based on the knowledge of the equation of state and on the requirement that the nuclear 
matter density is doubled for a total overlap of the two reacting nuclei. 



B. Folding Model Potential 

We consider two nuclei with one-body deformed densities pi and subjected to vibra- 
tional fluctuations, and center of masses separated by the distance r. Then the nuclear 
potential between these two nuclei can be evaluated as the double folding integral, 

yN{.r) = j dri j dr2 Pi(ri) p2(»'2) ^^(^'12), (11) 

where ri2 = r + r2 — ti. The central part of the effective NN interaction f (^12) contains 
no spin or spin-isospin terms since we analyze fusion reactions of spin nuclei, where the 
spin terms are relatively unimportant. We are then left with a direct part, 

Vdir[ri2) = ^^oo(?'i2) H 1 voi{ri2), (12) 

Ai A2 

which depends on isospin since N ^ Z in all the cases of interest. The exchange part, which 
takes into account the effect of antisymmetrization under the exchange of nucleons between 
the two nuclei, is modeled by the contact interaction, 

^ ^ ff , N,-Z, N2- Z2 J \ . 
VeA^u) = i^Joo + — j-^ j-^ — -^01 j S{ri2). (13) 

The density independent part of the effective nucleon-nucleon force that we use is the 



Reid parametrization of Michigan-3- Yukawa (M3Y) interaction 4S |. The explicit form of the 
expressions for t>oo, foi, -^00 and Jqi are given, for example, in Ref. 5(|- Moreover, we neglect 
the possible energy dependence of these parameters. For example, a variation 6E ^ 25 
MeV between the minimum and maximum bombarding energies considered in the fusion of 
^^Ni+^'^Ni gives a variation of 5Jqq ~ 0.27 MeV ■ fm^ MeV in the Reid parametrization, a 
values that is small compared to Jqo ~ —276 MeV ■ fm^. 

As we noted earlier, a potential that is based on the M3Y interaction predicts correctly 
the ion-ion potential for peripheral collisions. However, reactions that are sensitive to the 

10 



Nucleus 


rop (fm) 


Op (fm) 


ron (fm) 


an (fm) 


58Ni 


1.107 


0.4673 


1.0836 


0.5124 


64Ni 


1.0652 


0.575 


1.0852 


0.532 




1.1025 


0.535 


1.105 


0.575 



TABLE I: Parameters of the proton and neutron density distributions, Eq. p5|) . which have been 
used in the calculation of the direct and exchange, double-folding potential. 



potential at smaller distances are not reproduced [37|. To cure this deficiency, the ion- 
ion potential is supplemented with a short-ranged repulsive potential which, according to 
the discussion in the previous subsection, is proportional to the overlapping volume of the 
reacting nuclei. This is simulated in Eq. (fTTl) by using the effective contact interaction. 



V,ep{ri2) = Kep 5(ri2). (14) 

48| and use a relatively sharp density pro- 



We follow the procedure proposed in Ref. 
file, characterized by the diffuseness Orep, when calculating the repulsive potential from the 
double-folding integral, Eq. (|TT1) . The strength Vrep of the repulsive interaction is then ob- 
tained from the condition Eq. (fTOj) . which relates the total nuclear potential at the coordinate 
origin to the nuclear incompressibility. 



C. Parametrization of densities 

The densities that appear in Eq. (ITTl) are parametrized with Fermi-Dirac distribution 
functions 



Po 



(15) 



1 + exp[(r - i?oi)/a]' 

1/3 

where i?oi = ^oi^i ■ The parameters of the proton and neutron density distributions we 
have used are listed in Table [H The parameters for the proton density of ^^Ni are taken 
from a compilation of elastic electron scattering data [5l|. For the neutron distribution 
we choose the parameters such that they are in the range of what one would expect for a 
moderately neutron-rich nucleus, and such that the barrier of the M3Y+repulsion potential 
is close to the one predicted by the Akyiiz-Winther potential, since this potential gives a 
good description of the data in the barrier and high-energy region {4]. The parameters 
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Reaction 


Orep (fm) 


Frep (MeV-fm3) 


K (MeV) 


^pocket (MeV) 




0.341 


510.1 


234 


89.31 




0.405 


511.0 


228 


84.98 


64Ni+100Mo 


0.375 


488.7 


226 


117.8 



TABLE II: Parameters used in the calculation of the repulsive, double-folding potential for three 
fusion reactions, and the associated incompressibility and pocket energy. 




r (fm) 



FIG. 1: (Color online). Coulomb plus the nuclear M3Y+repulsion, double- folding potential for 
^^Ni-|-^^Ni. Results are shown for different values of the nuclear incompressibility K. 

for the proton and neutron densities in ^^Ni and ^°'^Mo we have used are identical or close 
to those predicted by the Hartree-Fock-Bogoliubov (HFB) method using the BSk2 Skyrme 



force 



52|. 



D. Heavy-ion potential study cases 



Applying the formalism described in subsections A and B, we vary the parameters Orep 
and Vrep of the repulsive potential so that Eq. (|TOl) is approximately fulfilled for the nuclear 
incompressibility K predicted by the Thomas- Fermi model 46||. The parameters we obtain 
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r (fm) 



FIG. 2: (Color online). Comparison of the double folding potentials with (solid curve) and with- 
out (long-dashed curve) repulsion, and the Akyiiz-Winther potential (short dashed curve) for the 
system ^^Ni+^^Ni. 

for the three heavy-ion systems we have studied are listed in Table [TTl together with the 
incompressibility of the compound nuclei obtained in the Thomas-Fermi model, and the 
energy of the resulting pocket that appears in the Coulomb plus nuclear potential. 

To illustrate the dependence of the total ion-ion potential on the nuclear equation of state 
we plot in Fig. [T]the Coulomb plus nuclear potential for the dinuclear system ^^Ni+^^Ni for 
different choices of the nuclear incompressibility K. The solid black curve corresponds to 
the incompressibility K ^ 234 MeV as inferred from the Thomas-Fermi model for ^^^Ba. 
The effect of K on the inner part of the barrier and the depth of the pocket is essential. A 
soft EOS provides a deep pocket whereas a hard one results in a shallow pocket. 

The effect of the repulsive core is displayed in Fig. [2j We compare the total poten- 
tials obtained from standard M3Y heavy-ion potential, the M3Y+repulsion potential for 
^^Ni+^^Ni, and the the Akyiiz-Winther, for the system ^^Ni+^^Ni. Due to the lower barrier 
and an abrupt decrease in the inner region, the M3Y potential cannot accurately reproduce 
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39| 



the data, a fact that already has been pointed out in the literature 

In Fig. [3] we show the profiles of the Coulomb plus nuclear potentials for four heavy-ion 
systems. The solid curves are based on the M3Y+repulsion potentials we use in this work, 
whereas the short-dashed curves are based on the AW potential. It is seen that the pockets 
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FIG. 3: (Color online). Ion-ion potentials for ^«Ni+^«Ni, '^''Ni+^^Ni, ^^Ni+^^'Ge, ^''Ni+^'^'^Mo. The 
solid curves are the potentials employed in the present work. The short-dashed curves are based 
on the Akyiiz-Winther potential, which was used in Refs. . The dashed strips show 

for each system the experimental boundaries of the threshold energy Es 



predicted for the MSY+repulsion potentials are much shallower than predicted by the AW 
potential, in particular for the lighter systems. The dashed region in Fig. [3] shows the energy 
Eg where the experimental S-factor develops a maximum, and the width of it illustrates the 
uncertainty in the extracted value of Eg jo]. 



V. COUPLED-CHANNELS APPROACH 



We use the same approach as in previous publications (see 531] and references therein), 
i.e. coupled-channels calculations performed in the so-called iso-centrifugal or rotating-frame 
approximation, where it is assumed that the orbital angular momentum L for the relative 
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motion of the dinuclear system is conserved. The rotating frame approximation allows a 
drastic reduction of the number of channels used in the calculations. 

The only ingredient we are going to change in the formalism is the ion-ion potential. 
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Like in the case of other nuclear potentials we resort to the "proximity" approximation 
which states that the heavy-ion potential is a function of the shortest distance between the 
nuclear surfaces of the reacting nuclei. Since for spherical nuclei, the relative orientation is 
not coming into play, the fragments have the natural tendency to be polarized along the 
collision axis, which is defined by the radial separation vector r. Surface vibrations are then 
taking place parallel (along) the direction r. In this case the shortest distance between the 
two surfaces is 



Ri-R2~6R, (16) 



where 



and r specifies the spatial orientation of the dinuclear system in the laboratory system and 
a^ji are the deformation parameters. In the rotating frame approximation, which we will 
use, the direction of r defines the z-axis. The only vibrational excitations that can take 
place are therefore the /i = components, since Y\^{z) oc 5^j,o- 

The cases under study in this paper refer only to the generation of vibrational excitations 
of surface modes. In previous papers it was shown that the inclusion of linear and quadratic 
interactions is necessary and often sufficient to fit the data at least in the intermediate 

and references therein). It is then necessary to consider only the 



energy region [see 



spherical part of the double-folding potential, and its first and second order derivatives. The 
nuclear potential for the elastic channel is then given by Eq. ( fTTl) . which is most conveniently 
calculated by the Fourier transform 54l |. 

V{.r) = ^ J dq pi{q) p2{q) v{q) Jo{qr). (18) 

Here Pi{q) is the Fourier transform of the spherical ground state density of ion i, v{q) the 
Fourier transform of the effective A^A^ interactions, and jo{qr) stands for the spherical Bessel 
function. This is actually the potential we have already shown in Figs. 1-3. 

The non-spherical part of the nuclear potential that results from the difference between 
the total interaction and the potential in the elastic channel is expanded up to second-order 
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in the surface distortion ([T 



(19) 



It is seen that the ground state expectation of this interaction, (0|(5VAr|0), is zero, but the 
second order term will give a non-zero contribution to the diagonal matrix element in an 
excited state. One can show that this prescription is exact for a harmonic oscillator up to 
second order in the deformation amplitudes. We include a similar expansion of the Coulomb 
field, SVc, but only to first order in the deformation amplitudes 53|. These expressions are 
inserted into the C. C. formalism in the rotating frame approximation which singles out only 
axially symmetric distortions (aA^=o), i-e.. 



2Mn 



^ L{L + 1) 



dr 



+ 



ni,n2 / 



{nin2 I 6Vc + 6Vn \ mim2)um^rn2{r), (20) 
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where E is the relative energy in the center of mass frame, L is the conserved orbital angular 
momentum, and Mq is the reduced mass of the dinuclear system. The C. C. equations f l20|) 
are written for two coupled vibrators of eigenenergy erii,n2 ^ind consequently the radial wave 
function u{r) is labeled by the quantum numbers ni and n2- Expressions for the matrix 
elements of 6V in the double-oscillator basis are given in jl5 |. 

The above C. C. equations are solved with the usual boundary conditions at large dis- 
tances and appropriate in-going-wave boundary conditions imposed inside the barrier, more 
precisely, at the radial separation where the potential pocket attains its minimum. The fusion 
cross section results then from the total in-going flux. The calculations include one-phonon 
excitations of the lowest 2^ and 3~ states in target and projectile, and all two-phonon and 
mutual excitations of these states up to a 7.2 MeV excitation energy. This energy cutoff was 
chosen so that all of the two-phonon states were included in the calculations for ^'^Ni+^^Ni 
and ^^Ni+^'^°Mo, whereas the two-phonon octupole states were excluded in the calculations 
for ^'^Ni+^^Ni. The necessary structure input for ^"^Ni and ^°''Mo is given in Ref. 55|; the 



input for Ni is from Ref. 



15|. 
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VI. CROSS-SECTIONS, 5-FACTORS AND LOGARITHMIC DERIVATIVES 



In what follows we analyze the fusion data for three systems for which there are strong 
indications of a hindrance, -or steep fallofF of the cross sections, at energies well below the 
Coulomb barrier. For completeness we also show in section VI. 4 the high-energy behavior 
of the cross sections. All calculations are based on the coupled-channels formalism outlined 
in the previous section, and further interpretations of the results are presented in Sect. VII. 



1. ^^Ni+^^m 

The first case we consider is the fusion reaction of ^^Ni+^^Ni. For this system data are 
available from an older experiment and the smallest cross sections are in the range of mb 
[3]. The data are compared in Fig. IH to several C. C. calculations that are based on the 
same structure input as in Ref. 15(1 and on the M3Y+repulsion potential, but they differ 
in the assumed value of the nuclear incompressibility. We see that for a hard nuclear EOS 
{K=290 MeV) the calculated low-energy cross sections deviate strongly from the measured 
values. Obviously in this case the origin of the mismatch is primarily due to the existence 
of a very shallow pocket which reduces the classically allowed region and hinders the fusion 
drastically. For a soft nuclear EOS (^=170 MeV) the potential pocket is lowered and the 
calculated cross sections approach the results calculated with the Akyiiz-Winther potential. 
The best fit to the data is obtained for the nuclear incompressibility 7^^=234 MeV, which is 
the value one obtains in the extended Thomas-Fermi model. 

Some of the results shown Fig. H] are repeated in Fig. [5], where they are compared to 
similar results for the fusion of ^"^Ni+^^Ni and ^^Ni-|-^°°Mo, which will be discussed in more 
detail below. The C. C. calculations for the M3Y+repulsion potential (solid curves) were 
all performed with the nuclear incompressibility derived from the extended Thomas-Fermi 
model. The long-dashed curves are the no-coupling limits (NOG) obtained with the AW 
potential. 

In Fig. [6] (left panel) we compare the experimental S'-factors for ^^Ni-|-^^Ni to the those 
calculated with the M3Y-|-repulsion and the AW potentials. In the last case the mismatch 
with the data is already evident starting at energies 4-5 MeV below the barrier. The former 
case, on the other hand, provides a good description of the available data and, most im- 
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FIG. 4: (Color online). Experimental fusion excitation function for the system ^^Ni+^^Ni [7[ is 
compared to various C. C. calculations that are based on the M3Y+repulsion potential (assuming 
different values of the nuclear incompressibility K), and on the AW potential. 

portantly, reproduces the trend of the experimental S'-factor to produce a maximum. The 
calculation also predicts a second maximum at even lower energies where no data exist. It 
remains to be seen whether future measurements will confirm the two predicted maxima. 

In Fig. [7] (lower panel) we plot the logarithmic derivative of the energy-weighted cross 
section for the case ^^Ni+^^Ni. We obtain a nice description of the available experimental 
data when we use the M3Y+repulsion potential (solid curve). As expected from the behavior 
of the S-factor, L{E) produces a maximum (just above the constant S'-factor limit) followed 
by a local minimum, before it finally diverges at lower energies. To ascertain this behavior, 
and also the predicted excitation function and S'-factor for ^^Ni-|-^^Ni, high-precision fusion 
data are necessary at very low bombarding energies, 

2. ^'^Ni+^^m 

Next we consider the system ^^Ni+®^Ni. This case was already discussed in our previous 
publication [Ij] but for the sake of completeness we find it necessary to recall it and add more 
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or online). Experimental fusion excitation functions for the systems ^^Ni+^**Ni 



_58i 



li B, 



Ni+'''*Ni and ''^Ni+^'^^Mo [5] are compared to various C. C. calculations described in the 
text, and to the no-coupling (NOC) limit for the AW potential. 



features like, for example, the logarithmic derivative. 

The system ®^Ni+^^Ni was advocated to be a very good choice for precise measurements 
since there is no contamination from reactions with heavier isotopes in the target, or from 
lower-Z isobaric components in the beam, which can dominate the yield at extreme sub- 
barrier energies jj]. The ATLAS data for this case are in a good agreement with previous 



results (56| at energies above the Coulomb barrier but not around and below the barrier. 
However, the agreement is good around and below the barrier with an older experiment j^]- 
It should be also mentioned that the ATLAS data provide cross sections down to tenths of 
nanobarns js] compared to the 0.3 mb that was reached in the older experiment 3]- 

To obtain the incompressibility of i^' ~ 228 MeV for total overlap of the reacting nuclei, 
as predicted by the Thomas- Fermi model, we use a strength of = 511 MeV and a 
diffuseness of the repulsive part of the density distribution of arep ~ 0.4 fm for both protons 
and neutrons. 

In Fig. [5] (middle panel) we compare the experimental excitation function for the fusion 
reaction ^^Ni + ^^Ni — >-^^^Ba with the C. C. results obtained using the Akyiiz-Winther 

n 

potential like in Ref. [2| (dotted line) and with the M3Y+repulsion potential (solid line). 
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FIG. 6: (Color online). The experimental S'-factors for the systems ^^Ni+^^Ni p;] (left panel), 
^^Ni+^^Ni {4] (middle panel) and ^^Ni+^'^'^'^Mo [5] (right panel) are indicated by solid circles. They 
are compared to the coupled-channels calculations performed with the MSY+repulsion (solid curve) 
and Akyiiz-Winther (dashed curve) potentials. 



The same recipe for the C. C. calculations was used in both cases. The diffuseness of the 
Akyiiz-Winther potential we use is a = 0.676 fm. One should also recall that calculations 
using a modified Akyiiz-Winther potential, where the diffuseness inside the barrier was set 
to ttj = 5 fm, were not able to the date to explain the steep fall-off of the measured cross 
section . This is not surprising because the modified AW potential has a very deep pocket 



(see Fig. 4 of Ref. 53| for an interior diffuseness of = 10 fm.) 

We conclude, from the inspection of Fig. O that the agreement with data, when using the 
M3Y-frepulsion potential, is sensitively better than the one provided by the Akyiiz-Winther 
starting at 90 MeV, and not only for the 4 lowest experimental data points. The excitation 
function obtained with the M3Y+repulsion potential has the right slope, not only because 
the potential attains a higher-lying pocket but also because the curvature of the barrier is 
different, with a thicker barrier in the overlapping region. Thus the best P^r point is only 

/N = 0.86. This value is obtained by applying the energy shift AE = 0.16 MeV to the 
calculated excitation function. The best fit obtained with AW potential, on the other hand, 
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FIG. 7: (Color online). Logarithmic derivatives of the energy-weighted cross sections for ^^Ni+^^Ni 
(lower panel), ^^Ni+^^Ni (middle panel), and ^^Ni+^'^'^Mo (upper panel). Experimental results 
derived from the cross sections shown in Fig. [S] are compared to coupled-channels calculations 
performed with the M3Y-|-repulsion (solid curve) and the Akyiiz-Winther (dashed curve) potential. 
The top dashed curve in each panel is the prediction for a constant S-factor. 
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gives a /N = 10 and requires an energy shift of AE = 0.9 MeV. 

The S-factor representation of the ^^Ni+^^Ni fusion cross sections is displayed in the 
middle panel of Fig. O The clear maximum in the measured S'-factor is reproduced only 
by the M3Y+repulsion potential. At this point we recall the experience gained in the 



past on molecular resonances. As shown in Ref. [23|], the S'-factor exhibits a sequence 
of quasi-molecular resonances, sandwiched between a limiting interior threshold and the 
Coulomb barrier. In the present case we obtain a maximum that is too broad to be assigned 
to a resonance, the curvature in the S'-factor being explained by the shallow pocket in 
the potential. The maximum of the theoretical curve corresponds approximately to the 
maximum of the experimental data. 

The logarithmic derivative of the energy-weighted cross section is shown in the middle 
panel of Fig. [71 and a nice agreement with the data is obtained in the C. C. calculations 
that are based on the M3Y-|-repulsion potential. 

3. ^-^Ni+^^^Mo 

For the repulsive part of the potential we choose a^ep = 0.375 fm and VJ-ep = 488.7 MeV, 
values that are matching the incompressibility K = 226 MeV of the compound nucleus 
^^^Yb. To improve the fit to the data we included up to three phonon excitations of the 
quadrupole mode in ^°°Mo using the structure parameters given in Ref. [55j]. However, the 
agreement with the data seen in the right panel of Fig. [5] is clearly not as good as in the 
other two cases shown in the same figure. The reason is that the C. C. effects are very 
strong for this heavy-ion system and the calculations have not fully converged with respect 



to multiphonon excitations, as discussed in Ref. [55|]. Another problem is that the nuclear 
structure properties of multiphonon states are often poorly known, so we will not try to 
improve the fit to the data here. 

For the S-factor we reproduce roughly, as can be concluded from an inspection of the 
right panel in Fig. [6l the trend to develop a well pronounced maximum at the lowest mea- 
sured energies and provide a reasonable estimation for this maximum. Also the logarithmic 
derivative manifests the tendency to develop a divergent behavior at low energies (Fig. [7l 
upper panel). The appearance of a local maximum at 124 MeV, on the other hand, is most 
likely caused by a poor convergence with respect to multiphonon excitations. 
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FIG. 8: (Color online). Linear plot of the fusion cross sections shown in Fig. [5] for the two systems 
s^Ni+^^Ni [4,] and ^'^Ni+^^Ni [T]. 

J^.. High-energy behavior of fusion cross sections 

The cross sections we obtain in coupled-channels calculations are suppressed when com- 
pared to the no-coupling limit. This can be seen in Fig. [5] but the suppression looks small 
on a logarithmic plot. We therefore show a linear plot in Fig. [8] of the same cross sections 
for the two nickel isotopes. In both cases we see that the coupled-channels calculations 
that are based on the AW potential are suppressed at high energies when compared to the 
no-coupling limit. This is a well-known C. C. effect but the suppression is not always large 
enough to explain the measurements, in particular for heavy systems lOj. 

We also see in Fig. 8 that the C. C. results we obtain with the M3Y-f-repulsion potential 
are suppressed even further, and the suppression increases with increasing energy. This is 
naively what one would expect with shallower pockets at higher angular momenta. It is 
unfortunate that the data considered here do not reach very large cross sections, so it is 
difficult to assess the significance of the suppression we obtain. However, the trend is very 
encouraging because it helps explaining why the high-energy fusion data considered in Ref. 
[lo| are suppressed so much compared to C. C. calculations and why it was necessary to use 
a large diffuseness of the ion-ion potential when fitting the data in the no-coupling limit. 
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VII. INTERPRETATION OF THE RESULTS 



In the previous section we obtained an excellent fit to the fusion data for the two nickel 
isotopes in the C. C. calculations that were based on the M3Y+repulsion potential and 
in-going-wave boundary conditions. Here we discuss the significance of these results and the 
implications for other observables. We have in this connections also tested other types of 



repulsive cores, like the gaussian one of the type used in [32] but we arrived at the conclusion 
that they are inadequate in reproducing the fusion data we have investigated here. Thus 
there appears to be some limitations or constraints in constructing the shape of the ion-ion 
potential inside the barrier. 

A. Models of fusion 

The description of heavy-ion fusion is often based on the assumption that fusion occurs 
as soon as the Coulomb barrier has been penetrated. More precisely, it is expected that 
the fusion cross section obtained from the absorption in a short-range imaginary potential 
is essentially identical to the cross section obtained with in-going-wave boundary conditions 
(IWBC) [57]. Moreover, the calculated fusion cross section is expected to be insensitive to 
small variations in the radius where the absorption or IWBC are imposed. 

In coupled-channels calculations, heavy-ion fusion is usually simulated by IWBC that 
are imposed at the minimum of the potential pocket, as for example in the computer code 
CCFULL 581] and also in the present work. The assumptions mentioned above have therefore 



rarely been challenged, and in cases where they have been tested, at energies close to the 



59|, 



Coulomb barrier, they have usually passed the test. One exception is discussed in Ref. 
where it was noticed that the calculated fusion cross section in one particular case showed 
an unacceptable variation with respect to the radius where the IWBC were iniposed. This 



occurred in a case with particularly strong couplings. Another example is Ref. 55| , where it 
was pointed out that certain unwanted oscillations in the cross sections, obtained at extreme 
subbarrier energies using the IWBC, could be suppressed by including a weak, short-range 
imaginary potential in the calculations. 

In this work we try to demonstrate that the fusion at extreme subbarrier energies is sensi- 
tive not only to the thickness and the height of the Coulomb barrier but also to the minimum 
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FIG. 9: (Color online). Dependence of the calculated fusion cross section (dashed curves) on the 
diffusenes aw of the imaginary potential. The results are compared to the data for ^''Ni+^^Ni J] 
and the calculation, which is based on the IWBC without any absorption (the solid curve). The 
triangular data point at the lowest energy is an upper limit. 



of potential pocket, where the IWBC are imposed. It is noted that the calculated fusion 
cross section, which is based on the IWBC, will vanish sharply as the energy approaches the 
minimum of the pocket. This sharp behavior is apparently what is needed to fit the data at 
extreme subbarrier energies, at least in the case of ^^Ni+^^Ni. 

We have also calculated the fusion of ^"^Ni+^^Ni using, in addition to the M3Y+repulsion 
potential, a short-range imaginary potential of the form proposed in Eq. (3) of Ref. 55|. 
The results for different values of the diffuseness aw are shown in Fig. [9] by the dashed 
curves. They all make an excellent fit to the data above 10 fih but they deviate significantly 
from the data at smaller cross sections. In contrast, the calculation discussed ealier (solid 
curve), which is based on the IWBC without any imaginary potential, is in perfect agreement 
with the data. 

There is a simple reason why the dashed curves are enhanced compared to the solid curve 
(IWBC) in Fig. [9] at the lowest energies. First of all, the fusion obtained from the IWBC 
without an imaginary potential can only occur in the elastic channel at low energy, say below 
87-88 MeV in the case considered here, and the cross section goes sharply to zero at the 
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minimum of the pocket which is 85.4 MeV. The inelastic channels are energetically closed to 
fusion at these low energies and the associated wave functions are (exponentially) decaying 
at the boundary. When an imaginary potential is turned on, the fusion/ absorption becomes 
possible through all channels. This explains qualitatively why the dashed curves are above 
the solid curve at low energy and why the conventional assumption about fusion described 
in the beginning of this subsection is not correct at extreme subbarrier energies. 



B. Average spin for fusion 

Another signature of a shallow pocket in the total ion-ion potential is a narrowing of the 
spin distribution for fusion as the center-of-mass energy decreases and approaches the pocket 
energy. In the past it has always been believed that the average angular momentum for fusion 
would approach a constant at low energy. This is the behavior that has been predicted by 



model calculations, including the C. C. [53|, but it has not really been tested experimentally. 
An example is shown in Fig. [101 where the measurements of the 7-ray multiplicity from the 
compound nucleus formed in the fusion of ^^Ni+^^Ni have been converted into an average 
angular momentum for fusion [56]. The thin dashed curve shows the prediction based on 
the AW potential in the no-coupling limit (NOC (AW)). It approaches a constant value at 
low energy but the data are always above that limit. The solid curve in Fig. [TO] shows the 
results we obtain in the C. C. calculations we discussed earlier, which were based on the 
MSY-f-repulsion potential. It is seen that these calculations predict a narrowing of the spin 
distribution as the energy approaches the pocket energy. 

After completing our studies we realized that the low-energy behavior of several observ- 
ables can have a strong sensitivity to the couplings to multi-phonon states. Although the 
fusion only occurs in the elastic channel at energies close to the minimum of the pocket, 
the polarization of inelastic channels can still have a large effect. We found, in particular, 
that couplings to the two-phonon octupole states are very important. This is illustrated in 
Fig. [To] by the dotted curve which was obtained without any couplings to the two-phonon 
octupole states. It is seen that the calculation in this case develops a rather sharp peak 
at 87.7 MeV. The peak disappears when the coupling to the two-phonon octupole states is 
included, as illustrated by the solid curve. Unfortunately, the data cannot tell us which of 
these two calculations is the most realistic. 
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FIG. 10: (Color online). Average angular momentum for the fusion of ^^Ni+^^Ni obtained in 
C. C. calculations based on the M3Y+repulsion potential with (solid curve) and without (thick 
dotted curve) the effect of couplings to the two-phonon octupole states. The thin dashed curve 



was obtained in the no-coupling limit (NOC) using the AW potential. The data are from 



56| 



C. Structures in calculated cross sections 



The local maximum that appears in the thick dotted curve in Fig. [TO] at low energy 
is not a resonance. It is the result of a vanishing wave function in the elastic channel 
at the radial separation where the IWBC are imposed. When this condition is fulfilled 
at low energy, where fusion is restricted to the elastic channel, the fusion probability will 
vanish. When it occurs for a range of low angular momenta, it will result in a large average 
angular momentum for fusion and that explains the appearance of the peak in Fig. [TOl 
The suppression of the fusion probability at low energy (which occurs when the two-phonon 
octupole states are not included in the calculation) can also produce a local minimum in the 
^-factor. This is illustrated by the dotted curve in Fig. [TTl The solid curve, which includes 
the effect of couplings to the two-phonon octupole states, shows a single, broad maximum. 
In this particular case there is a clear preference for the solid curve which makes an excellent 
fit to the data. 

In the calculations we presented for ^^Ni+^^Ni we did not include the two-phonon octupole 



27 



10^ 



> 

(U 



I 



"a, 

X 

p-l 



10^ 



10' 



10" 



CC (M3Y+rep, no 2PH(3")) 
CC (M3Y+rep) 





84 86 88 90 92 94 

E (MeV) 



FIG. 11: (Color online). S factors for the fusion of ^''Ni+^^Ni obtained in the C. C. calculations 
based on the M3Y+repulsion potential with (solid curve) and without (dotted curve) the effect of 
couplings to the two-phonon octupole states. The data are from Ref. 14]; the triangle is an upper 
limit. 

states in the C. C. calculations because the excitation energy is very large, almost 9 MeV, 
so we did not expect these states would be important. In view of the above discussion 
it is now understandable why the calculated S-factor for ^^Ni+^^Ni shown in left panel of 
Fig. develops a double-peaked structure at low energy. We have therefore repeated the 
calculations and included the two-phonon octupole states. The resulting S'-factor exhibits a 
single, broad peak, just as we saw in the fusion of ^^Ni+^^Ni. 

It is not clear why the couplings to the two-phonon octupole states play such an impor- 
tant role as discussed above. However, the analysis of the ^^Ni+^^Ni data shows a perfect 
agreement with a single, broad S factor peak, whereas the analysis of the existing fusion 
data for ^^Ni+^^Ni fusion has a strong preference for the double-peaked S'-factor curve. In 
view of these findings it is of interest to continue the fusion measurements for ^^Ni+^^Ni to 
even lower energies because the existing data [3] shown in the left panel of Fig. [6] do not 
verify explicitly the double-peaked structure. 
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FIG. 12: (Color online). Plot of Vpocket (triangles) and Eg (solid points) versus ZxZ^^fJi. for the 
three reactions : ^^Ni+^^Ni, ^^Ni+^^Ni, ^^ni+ioomo. 



D. Correlation between the pocket energy and Eg 



As we mentioned in the Introduction, Ref. 



3| reported a correlation between the energy 
where the S'-factor has a maximum and the parameter ZiZ2y/JI, with /i = AiA2/{Ai + A2). 
In view of the idea of a shallow pocket advocated in the present paper we considered it 
useful to study the dependence of the minimum of the potential Ipocket, along with Eg, on 
the parameter ZiZ2^fJi. This is illustrated in Fig.[12]for the three cases we have investigated 
in this work. We note the pocket energy V^ocket follows closely the energy Eg, where the S 
has a maximum, the two results being separated by about 3 to 5 MeV. 



VIII. CONCLUSIONS 



We conclude that in order to explain the experimental fusion data at extreme subbar- 
rier energies it is not necessary to adopt the hypothesis of an abnormal diffuseness of the 
internuclear potential as advocated in some publications. Such an assumption would be 
equivalent, within the folding model, to an unusual diffuseness of the target and projec- 
tile nuclear matter distributions. Proton diffusivities, as inferred from electron scattering 
measurements (^^Ni) or from HFB calculations (^^Ni and ^°°Mo) are able to reproduce the 
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excitation functions. Using a diffuseness of the neutron matter distribution in the range of 
0.5-0.6 fm, which is a reasonable assumption for the moderately neutron rich nuclei con- 
sidered in this paper, is necessary not only to realize a shallow pocket but also to obtain a 
good value for the height of the spherical Coulomb barrier. With these remarks we do not 
want to discard the appearance of a large diffuseness during the fusion process. Physically 
we expect an increase in the nuclear skin to appear in the neck (overlap) region but not over 
the entire surface of the reacting nuclei. 

We have tried to advocate in this work that the understanding of the experimental data 
requires more than simply a modification in the curvature of the barrier (as used in the 
Hill-Whcclcr approximation of the Wong formula), or the introduction of a simple recipe for 
the repulsive potential like a gaussian. Wc used the double-folding potential that is based on 
the Reid parametrization of the M3Y interaction, and realistic parameters of the proton and 
neutron distributions of both target and projectile. We have supplemented this potential 
with a repulsive potential that takes into account the incompressibility of the nuclear matter. 

We also conclude that it is necessary to define the fusion in terms of in-going-wave 
boundary conditions. When these consitions are imposed at the minimum of the potential 
pocket, it is possible to reproduce the steep energy dependence of the fusion data at extreme 
subbarrier energies. Simulating the fusion by the absorption in an imaginary potential, on 
the other hand, does not allow us to reproduce the steep falloff of data at the lowest energies. 

In the two cases that recently were measured at ATLAS, namely, ^^Ni-|-^^Ni and 
^^Ni-|-^°°Mo, our calculations show a single broad maximum in the -S-factor at the low- 
est measured energies in agreement with the measurements. In contrast, our calculations 
for ^^Ni+^^Ni show a doublc-pcakcd structure of the S-factor. Wc find that the predicted 
shape of the low-energy ^'-factor is very sensitive to the couplings to the two-phonon octupole 
states. Thus if we ignore these couplings, the S'-factor develops a double-peaked structure at 
low energy, whereas a strong coupling to the these states (as in a harmonic vibration) tends 
to produce a single, broad low-energy peak. It is therefore important to test the predicted 
S'-factor experimentally at lower energies, in particular in the case of ^^Ni-|-^^Ni. 

We noticed for all three fusion reactions studied in this paper a nice correlation between 
the minimum of the potential pocket Knin and the experimentally extracted reference energy 
Eg, where the S'-factor reaches a maximum. A systematic study of fusion reactions over 
a wider range of values of the parameter Z\Z2^lJi is necessary to confirm this apparent 
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correlation, and the conjecture of a repulsive core for overlapping configurations. We also 
mentioned another possible signature of a shallow pocket, which is the narrowing of the spin 
distribution at energies below Eg. This is what our C. C. calculations predict and it can be 
tested by measurements of the 7-ray multiplicity emitted from the compound nucleus. 

Before submitting this paper further evidence on the hindrance in the sub-barrier fusion 
of ^^Ca+^^Zr was reported 60|- This conclusion resulted from the analysis of the logarithmic 
derivative which exhibits a steep increase at the lowest measured cross sections. 
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